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ABSTRACT 


When  an  oscillating  piston  forces  the  gas  in  a  closed 
tube  to  vibrate  at  a  finite  amplitude  near  an  acoustic  reso¬ 
nance  frequency,  shock  waves  are  generated  which  travel  back 
and  forth  in  the  tube.  The  behavior  of  these  waves  in  tran¬ 
sition  from  the  acoustic  into  the  shock  region  for  air  at  stan¬ 
dard  temperature  and  pressure  was  studied  at  fundamental  reso¬ 
nance  frequencies  up  to  1000  cycles.  An  electrodynamic  vibra¬ 
tor,  capable  of  variable  acceleration  up  to  140  meters  per 
second  per  second,  was  used  as  a  piston  driver;  it  produced 
acoustic  pressures  in  the  tube  up  to  1 60  db  sound  pressure 
level . 

The  dependence  of  pressure  wave  shape  on  piston  motion 
was  investigated;  it  was  discovered  that  over  a  large  range  of 
piston  accelerations,  nearly  identical  wave  shapes  were  pro¬ 
duced  at  all  frequencies  provided  that  the  same  piston  accele¬ 
ration  was  maintained. 

Although  the  largest  shocks  occur  when  the  frequency  is 
adjusted  so  that  the  shocks  arrive  at  the  piston  just  as  the 
piston  reaches  its  maximum  withdrawn  position,  shocks  can  be 
observed  at  appreciable  fractions  of  a  period  earlier  or  later. 
The  effect  of  small  changes  in  frequency  away  from  resonance 
on  this  phase  difference  was  studied,  and  a  similarity  para¬ 
meter  was  determined  for  this  variation. 

Shock  strength  was  measured  as  a  function  of  both  the 
phase  of  shock  occurrence  at  the  piston  and  piston  accelera¬ 
tion.  It  was  also  found  that  acoustic  resonance  and  shock 
resonance  for  a  given  tube  occur  at  exactly  the  same  frequency. 
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1 .  Introduction 

When  an  oscillating  piston  forces  the  gas  in  a  closed 
tube  into  vibration,  a  number  of  interesting  effects  occur. 

If  one  undertakes  to  study  the  motion  of  the  gas  in  such  a 
tube  according  to  acoustic  theory,  results  are  obtained  which, 
although  reasonably  correct  in  most  instances,  predict  particle 
velocities  and  pressures  which  tend  to  infinity  at  resonance. 

Consider  a  cylindrical  tube,  closed  at  one  end  by  a  rigid 
wall  and  at  the  other  end  by  a  vibrating  piston.  Assume  that 
the  diameter  of  the  tube  is  sufficiently  small  so  that  the 
waves  travel  down  the  tube  with  plane  wave  fronts.  With  particle 
velocity,  u,  as  the  dependent  variable,  the  wave  equation  is 

d2u  _  j_  gfa 

9%2  c2  at2 

whose  solution  may  be  written 


where  cj  =  kc,  and  k  =  2Tt/a  ,  with  co  the  angular  frequency 
and  X  the  wavelength.  The  boundary  conditions  may  be  approxi¬ 
mated  by 

it 

At  x  =  0,  u(0,t)  =  uQe  ,  so  that 


U+  +  U_  =  U0 

at  x  =  L,  u(L,t)  =  0,  so  that 


—  ik  4  ikLi 

u+e  -4-  u_  e 


o 


(2) 


(3) 


solving 


u(x>t)  ^  Uc 


_  iuit 

e  .  (4) 


Sin  kj  L 

If  position  is  fixed  at  some  particular  value  of  x  and  if  L  is 
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held  constant,  then,  as  frequency  is  varied,  both  the  numerator 
and  denominator  of  (4)  will  vary.  When  kL  Is  some  multiple 
of  TT ,  the  particle  velocity  becomes  Infinite,  except  at  points 
where  k(L  -  x)  Is  also  a  multiple  of  TT  #  This  leads  to  an  ex¬ 
pression  for  the  length  of  a  closed  tube  at  resonance: 


n  A 
Z 


It  may  be  verified  by  experiment  that  even  though  the  ampli¬ 
tude  of  oscillation  of  the  driving  piston  Is  Infinitesimal 
compared  to  the  length  of  the  tube,  large-amplitude  velocity 
and  pressure  vibrations  are  obtained  when  the  length  of  the 
tube  Is  as  given  above.  It  Is  to  be  noted,  however,  that  the 
approximate  expression  given  by  (4)  becomes  Indeterminate  at 
the  ends  of  a  resonating  tube.  B.  Miller  and  L.  Olsen  [12] 
make  a  more  nearly  exact  analysis  based  on  the  wave  equation 
written  with  a  velocity-dependent  damping  term.  They  develop 
an  expression  for  the  length  of  a  resonant  tube  given  by 


where  c*  is  the  damping  constant.  In  the  usual  experimental 
situation  encountered  in  acoustics,  the  second  term  is  negli¬ 
gible  in  comparison  with  the  first,  and  the  expression  reduces 
to  that  given  in  the  approximate  analysis. 

In  practice,  of  course,  although  particle  velocities  and 
pressures  become  very  large  at  resonance,  the  effects  of  vis¬ 
cosity  and  heat  dissipation  are  such  as  to  limit  the  amplitudes 
to  finite  values.  If  the  piston  amplitude  is  truly  infinitesi¬ 
mal,  the  wave  shape  at  resonance  retains  its  sinusoidal  form, 
and  the  wave  may  be  referred  to  as  an  acoustic  wave. 
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As  the  amplitude  of  piston  oscillation  is  increased, 
the  behavior  of  the  gas  changes.  The  most  striking  pheno¬ 
menon  occurs  when  the  frequency  of  vibration  is  near  one  of 
the  acoustic  resonance  frequencies.  This  is  the  formation  of 
shock  waves  which  travel  back  and  forth  in  the  tube.  With 
these  finite  amplitude  vibrations,  the  pressure  and  velocity 
waves  make  a  radical  departure  from  the  sinusoidal  character 
observed  at  smaller  amplitudes,  and  it  becomes  obvious  that 
the  acoustic  theory  represented  by  (4)  no  longer  applies. 

Up  to  the  present  time,  a  complete  theory  describing 
such  finite  amplitude  vibrations  has  not  been  formulated. 
Several  investigations  of  varying  thoroughness  have  been  con¬ 
ducted,  which  have  led  to  somewhat  inconclusive  results,  and 
it  appears  that  more  experimental  data  are  needed.  Previous 
experiments  were  conducted  using  mechanically  driven  pistons 
which  of  necessity  moved  with  constant  amplitude,  since  with¬ 
out  special  linkages  it  is  somewhat  awkward  to  change  the  ampll 
tude  of  such  a  driver.  The  experiments  which  are  the  subject 
of  this  thesis  were  conducted  using  an  electrodynamic  vibrator 
which  was  easily  variable  over  a  large  range  of  amplitude  and 
frequency,  so  that  measurements  could  be  made  under  many  dif¬ 
ferent  conditions.  This  study  was  conducted  toward  a  threefold 
purpose: 

a.  To  show  the  effects  of  variation  of  amplitude  on  the 
occurrence  and  character  of  periodic  shock  waves. 

b.  To  increase  our  knowledge  of  the  effects  of  other 
parameters . 
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c.  To  demonstrate  the  feasibility  of  the  use  of  an 
electrodynamic  shaker  for  generating  shocks  for  further  ex¬ 
periment. 
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2.  Historical  Background 

Investigations  into  the  phenomena  of  nonlinear  gaseous 
vibrations  in  a  closed  tube  were  begun  in  1935  by  E.  Schmidt 
jl  ]  and  0.  Mayer-Schuchard  [2],  who  were  interested  in  deter¬ 
mining  the  nature  of  this  type  of  vibration  as  experienced  in 
Internal  combustion  engines.  They  used  a  tube  (length  1200  cm, 
radius  3.4-  cm,  piston  stroke  3.4  cm)  filled  with  air  at  nor¬ 
mal  density.  A  motorcycle  engine  was  driven  by  an  electric 
motor  and  its  piston  furnished  excitation  to  the  tube.  Near 
the  acoustic  resonance  frequencies,  pressure-time  histories 
measured  at  the  closed  end  of  the  tube  showed  marked  discon¬ 
tinuities.  Between  these  pressure  jumps  the  pressure  records 
followed  a  curve  resembling  a  sinusoid.  Schmidt  recognized  the 
jumps  as  shock  waves. 

Using  the  same  apparatus,  E.  Lettau  [3]  in  1939  made 
numerous  observations  of  both  pressure  and  velocity  at  several 
stations  in  the  tube.  His  measurements  were  made  by  electro¬ 
mechanical  means,  and  his  records  we re  obtained  by  photograph¬ 
ing  the  trace  made  by  a  moving  light  beam.  He  found  that  the 
shape  of  the  velocity  and  pressure  profiles  was  strongly  de¬ 
pendent  on  the  location  in  the  tube  at  which  the  measurement 
was  made  (see  Fig.  1),  and  postulated  an  explanation  involv¬ 
ing  the  summation  of  two  partial  waves  traveling  in  opposite 
directions. 

Attempts  to  solve  for  the  motion  of  the  gas  in  a  closed 
tube  by  mathematical  means  were  begun  somewhat  later.  J.  Kel¬ 
ler  [4],  in  1952,  obtained  exact  solutions  for  the  one-dimen- 
eional  gas  dynamic  equations,  assuming  an  inviscld,  non-heat- 
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PRESSURE 


PARTICLE 

VELOCITY 


PARTICLE 

DISPLACEMENT 


Figure  1 .  Variation  over  one  cycle  of  pressure, 

particle  velocity  and  particle  displace¬ 
ment  at  different  stations  in  a  closed 
tube  resonating  at  its  fundamental  fre¬ 
quency.  (After  Lettau  [3]  ) 


conducting  medium  with  a  particular  equation  of  state,  often 
referred  to  as  a  Chaplygin  gas.  One  of  the  results  obtained  by 
Keller  is  that  shocks  should  not  occur  in  the  tube  unless  the 
piston  amplitude  is  in  excess  of  A/2  ,  where  A  is  the  wave¬ 
length.  That  experiment  produces  shocks  at  piston  displace¬ 
ment  amplitudes  of  considerably  less  than  this  magnitude  shows 
that  in  this  respect,  at  least,  Keller's  solution  does  not 
accurately  describe  the  phenomenon. 

R.  Betchov  [6]  made  a  theoretical  analysis  of  the  same 
problem,  this  time  assuming  a  tube  filled  with  a  perfect  gas, 
and  taking  into  account  the  effects  of  viscosity  and  heat  con¬ 
duction.  He  chose  to  make  linear  approximations  of  the  motion 
of  the  gas,  which  he  stated  are  valid  in  every  portion  of  the 
tube  except  in  the  vicinity  of  the  shocks  themselves.  A  dis¬ 
continuity  between  the  linear  approximations  represents  the 
existence  of  shocks.  Prom  his  solution  he  prepared  charts  of 
the  variation  of  velocity  and  density  which  show  qualitative 
agreement  with  the  results  of  Lettau. 

R.  Saenger  and  G.  Hudson  [7],  [8  3  made  a  complete  theore¬ 
tical  analysis  based  on  the  idealized  physical  model  of  a  gas 
element  shown  in  Fig.  2.  This  shows  the  element  together  with 
the  boundary  layer  of  gas  between  it  and  the  tube  wall.  The 
element  moves  as  a  unit  under  the  action  of  pressure  forces  on 
its  opposite  faces  and  frictional  retarding  forces  at  the  rim. 
The  wall  acts  as  a  heat  reservoir  and  the  use  of  Newton's  law 
of  cooling  permits  the  development  of  one-dimensional  energy 
balance  relationships.  Saenger  and  Hudson  assumed  that  gas 
motion  at  resonance  may  be  represented  by  the  sum  of  a  finite 
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Figure  2.  Idealized  physical  model  of  gas  element 

used  by  Saenger  and  Hudson  for  their  theo¬ 
retical  development. 

continuous  component  and  a  finite  discontinuous  one,  both  of 
which  are  periodic  in  time.  These  were  obtained  using  the  non¬ 
linear  expressions  derived  from  the  basic  flow  equations,  the 
Rankine-Hugoniot  shock  conditions,  and  the  boundary  conditions 
on  the  medium.  They  were  then  used  to  make  a  theoretical  deter 
mination  of  shock  strength  and  mean  temperature  as  functions 
of  piston  amplitude,  tube  length  and  dissipation  parameters. 

Comparison  between  theory  and  the  meager  amount  of  experi¬ 
mental  data  available  at  the  time  showed  qualitative  agreement. 
One  difficulty  which  faced  Saenger  and  Hudson  was  the  determina 
tlon  of  the  proper  dissipation  parameters  for  heat  conduction 
and  frictional  effects.  They  were  unable  to  do  this  theoreti¬ 
cally,  and  were  forced  to  use  empirical  values  to  obtain  agree¬ 
ment  between  theory  and  experiment. 

Among  the  conclusions  obtained  in  their  treatment  is  the 
following:  for  piston  amplitudes  less  than  a  critical  value, 

periodic  shocks  cannot  exist.  Thus  for  a  tube  of  fixed  length 
driven  at  piston  amplitudes  below  this  critical  value,  it  may 
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be  expected  that  a  continuous,  though  not  necessarily  acoustic, 
treatment  is  appropriate.  This  implies  that  the  transition  from 
the  acoustic  region  to  the  shock  region  should  be  a  sudden  pro¬ 
cess,  with  the  shock  suddenly  appearing  fully  developed  in  place 
of  a  wave  of  acoustic  proportions. 
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3.  Description  of  Periodic  Shock  Wave  Phenomena 

In  order  to  observe  shock  wave  phenomena  in  a  closed  tube, 
it  is  necessary  first  to  match  the  driver  frequency  with  the 
length  of  the  tube  for  acoustic  resonance.  This  length  is  some 
Integral  multiple  of  A/2.  For  small  values  of  piston  dis¬ 
placement  only  the  pure  tone  characteristic  of  acoustic  vibra¬ 
tions  is  audible.  But  as  the  amplitude  is  increased,  a  buzzing 
sound  is  superimposed  which  indicates  the  presence  of  periodic 
shock  waves.  If  the  pressure-time  history  is  observed  on  an 
oscilloscope,  one  notices  the  steep  discontinuities  occurring 
at  the  same  time  in  each  cycle.  If  the  length  of  the  tube  is 
varied  by  a  few  millimeters,  or  the  frequency  by  a  few  tenths 
of  a  cycle,  the  shape  of  the  pressure  curve  changes  drasti¬ 
cally,  and  the  shock  wave  changes  in  magnitude.  If  either  fre¬ 
quency  or  length  are  changed  by  more  than  minimum  values,  the 
shock  condition  disappears  and  the  tube  returns  to  an  acoustic 
mode  of  vibration,  which  gradually  diminishes  in  strength  as 
the  tube  departs  further  from  the  resonant  condition. 

Fig.  3  is  an  x  -  t  diagram  showing  the  motion  of  the  shock 
wave  as  it  travels  in  the  tube,  with  respect  to  the  motion  of 
the  piston.  From  this  it  can  be  seen  that  the  shock  wave  reaches 
the  piston  at  the  instant  that  it  begins  the  lnstroke  portion 
of  its  cycle.  In  actuality,  the  frequency-tubelength  combina¬ 
tion  may  be  adjusted  so  that  the  arrival  of  the  shock  at  the 
piston  is  considerably  out  of  phase  with  respect  to  this  maxi¬ 
mum  outstroke  position.  It  must  be  mentioned  that  the  use  of 
straight  lines  on  the  diagram  of  Fig.  3  is  probably  an  ideali¬ 
zation,  since  it  is  likely  that  they  should  be  curved,  repre- 
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Figure  3.  Space-time  diagram  for  first  resonance  (f^  = 

c/2L)  showing  the  characteristic  lines  of  slope 
=  1/c  along  which  periodic  shocks'  travel. 

senting  changes  in  shock  velocity  with  local  conditions  [^.6 J  . 

Particle  velocity  measurements  can  also  be  used  to  show 
the  existence  of  shock  waves,  as  can  be  seen  from  Fig.  1  . 

The  apparatus  used  in  the  present  study  proved  to  be  very 
well  suited  to  the  study  of  periodic  shock  waves,  primarily 
because  it  allows  the  precise  setting  of  parameters  and  the 
quick  determination  of  their  effects  on  the  shock  waves.  The 
results  obtained  can  be  more  easily  understood  in  terms  of  the 
following  parameters  and  concepts. 

a.  Resonance 

The  concept  of  resonance  is  basically  an  acoustical  one, 
indicating  the  condition  at  which  the  tube  length  is  precisely 
adjusted  to  be  nA/2,  causing  the  sound  waves  in  the  tube  to  be 
of  maximum  strength.  How  the  acoustic  resonant  condition  is 
related  to  the  production  of  maximum  shock  strengths  had  not 
been  determined  prior  to  the  present  study. 
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b.  Piston  Motion 

In  the  process  of  making  theoretical  developments  concern¬ 
ing  periodic  shock  waves,  previous  writers  have  expressed  re¬ 
sults  in  terms  of  piston  amplitude.  Keller  and  Betchov  use 
piston  displacement  amplitude  in  its  dimensional  form,  while 
Saenger  prefers  to  use  a  non-dimensionalized  form  which  is  essen 
tially  proportional  to  piston  velocity  amplitude.  The  results 
of  Saenger  are  derived  in  terms  of  the  quantity  X/L,  where  X 
is  piston  displacement  amplitude  and  L  is  the  length  of  the 
tube,  seeking  to  make  theoretical  predictions  which  are  appli¬ 
cable  to  any  length  of  tube,  so  long  as  it  is  resonant  to  the 
driving  frequency.  But 

A= 

and  for  a  resonant  tube, 

i  =  2i  -  irc 

L  z  ~  u> 

Hence, 

X_  ___  OJX  _  _V_ 

T7  ttc  t rc  (5) 

showing  that  the  parameter  X/L  is  actually  proportional  to  pis¬ 
ton  velocity  amplitude.  Experimentally,  the  measurement  of  dis¬ 
placement  amplitude  presents  some  difficulty,  and  it  is  often 
easier  to  measure  piston  velocity  and  integrate  to  obtain  dis¬ 
placement  amplitude,  which  is  a  simple  matter  if  the  piston  mo¬ 
tion  is  sinusoidal. 

Considerable  evidence  was  obtained  from  the  present  study 
that  piston  acceleration  is  a  more  appropriate  parameter.  Ac¬ 
celeration  can  also  be  calculated  simply  from  piston  velocity 
measurements. 
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c.  The  parameter  A- 

One  of  the  most  critical  variables  affecting  the  strength 
of  the  shocks  and  the  shape  of  the  waves  is  the  phase  difference 
between  the  time  of  the  arrival  of  the  shock  at  the  piston  and 
the  time  that  the  piston  reaches  its  maximum  outstroke  position. 
This  parameter,  A,  is  illustrated  in  Fig.  4.  In  this  illustra¬ 
tion  the  smooth  curve  is  a  velocity  output  from  an  electromag¬ 
netic  pickup  attached  to  the  piston  armature.  When  this  curve 
crosses  the  horizontal  axis,  the  piston  is  at  one  of  its  ex¬ 
treme  positions.  The  phase  of  the  extreme  position  can  be  mea¬ 
sured  with  respect  to  the  occurrence  of  the  shock.  In  the  pre¬ 
sent  experiment,  the  measurement  of  velocity  was  made  at  the 
piston,  and  the  measurement  of  pressure  at  the  opposite  end  of 
the  tube.  One  must  assume  that  the  shock  takes  half  a  period 
to  travel  the  length  of  the  tube  in  each  direction.  The  para¬ 
meter  A  has  been  assigned  a  negative  sign  if  it  is  such  that 
the  piston  reaches  its  maximum  outstroke  position  before  the 
shock  appears  at  the  piston,  and  a  positive  sign  if  the  piston 
extremum  occurs  after  the  shock. 

d.  Tube  shape  and  dimension. 

It  seems  likely  that  tube  shape  and  dimension  would  have 
considerable  effect  on  the  behavior  of  the  resonating;  air  col¬ 
umn.  The  theory  of  Saenger  and  Hudson  involves  the  effects  of 
viscosity  and  heat  conduction,  both  of  which  would  be  sensitive 
to  changes  in  the  tube  dimension.  If  schlieren-optical  tech¬ 
niques  are  to  b#  used,  it  would  be  necessary  to  change  the 
shape  of  the  tube  to  a  rectangular  cross-section.  In  *the  pre¬ 
sent  study,  only  one  tube  was  used,  and  that  of  circular  shape. 
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4.  Description  of  Equipment 

a.  Piston  driver. 

For  the  study  of  threshold  behavior  of  periodic  shock  waves, 
the  desirable  characteristics  of  a  driver  include  (1)  suffi¬ 
ciently  high  amplitude  for  the  production  of  shocks,  (2)  con¬ 
tinuously  variable  frequency  over  a  wide  range,  (3)  easily  ad¬ 
justed  and  continuously  variable  piston  amplitude,  and  (4)  pure 
sinusoidal  motion.  Although  it  is  theoretically  possible  to 
obtain  all  of  these  qualities  in  a  mechanical  driver,  design 
difficulties  arise  which  would  ultimately  result  in  an  expen¬ 
sive  and  cumbersome  apparatus  which  would  be  difficult  to  use. 

The  electrodynamic  vibration  shake  table  was  decided  upon 
as  a  driver  closely  approaching  the  four  desirable  states  men¬ 
tioned  above.  For  this  study  the  equipment  used  was  a  Calidyne 
Model  A88  Shaker,  manufactured  by  the  Calidyne  Company,  Winches¬ 
ter,  Mass.  This  shaker  is  designed  with  the  objective  of  pro¬ 
viding  a  light  weight,  yet  rugged  armature  frame  for  transmit¬ 
ting  vibratory  forces  to  an  attached  load.  A  set  of  trunnions 
allows  it  to  be  tilted  for  horizontal  operation.  Tapped  holes 
are  provided  for  connecting  the  load  to  the  armature  of  the 
shaker.  Into  one  of  these  holes  a  piston  was  screwed  which  was 
used  to  excite  the  air  column  in  the  tube. 

Included  with  the  shake  table  were  a  Power  Supply  Model  68, 
an  audio  oscillator,  a  Power  Amplifier  Model  68,  and  a  D.  C. 

Field  Supply  Model  102C,  all  Calidyne  equipment.  FREQUENCY 
and  AMPLITUDE  controls  on  the  audio  oscillator,  as  well  as  a 
GAIN  control  on  the  power  amplifier,  allowed  variation  of  fre¬ 
quency  from  5  cps  to  10  kc,  and  of  rms  acceleration  from  zero 
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to  about  140  meters  per  second  per  second.  In  addition,  choice 
of  proper  impedance  match  (using  IMPEDANCE  MATCH  selector  on 
Amplifier  Model  68)  produced  sinusoidal  piston  motion  with  less 
than  one  percent  harmonic  distortion  in  the  velocity  signal  at 
frequencies  up  to  2000  cycles.  Specifications  of  the  Model  A88 
shaker  are  included  in  Appendix  I. 

The  Model  A88  Shaker  is  equipped  with  a  velocity  pickup 
which  consists  of  a  coil  attached  to  the  table  of  the  shaker 
so  as  to  move  in  the  field  of  a  permanent  magnet  attached  to 
the  shaker  body.  The  voltage  produced  is  proportional  to  the 
relative  velocity  between  the  shaker  table  and  body.  This  vol¬ 
tage  can  be  continuously  monitored  as  an  indication  of  the  ab¬ 
solute  table  velocity.  The  use  of  an  a .  c.  voltmeter  allows 
the  shaker  to  be  adjusted  for  operation  at  desired  velocity  am¬ 
plitudes. 

b.  Microphone 

An  Altec  21-BR-200  microphone  system  with  Amplifier  Type 
526b  was  used  as  a  pressure  pickup.  This  microphone  was  used 
because  of  its  rugged  construction  and  essentially  flat  frequen 
cy  response  curve.  Its  relatively  slow  rise  time  of  about  60 
microseconds  placed  a  limitation  upon  the  accuracy  of  observa¬ 
tions  made  at  higher  frequencies.  A  comparison  calibration  of 
the  microphone  pressure  sensitivity  was  performed  using  a  Wes¬ 
tern  Electric  microphone  Model  WE-640-AA  as  a  standard.  The  re 
suits  of  this  calibration  are  contained  in  Appendix  II. 

c.  Tube  and  Piston. 

The  study  was  performed  using  a  1-3/4  inch  diameter  alumi¬ 
num  tube  approximately  seven  feet  long.  A  simple  aluminum  pis- 
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ton  of  approximately  the  same  diameter  was  machined,  and  was 
connected  by  a  short  threaded  rod  to  the  shaker.  Initially,  a 
metal  bracket  was  constructed  Into  which  the  tube  was  to  be 
threaded;  the  bracket  was  to  be  bolted  to  the  face  of  the  shaker 
with  the  Intention  of  holding  the  end  of  the  tube  rigidly  In 
place.  Because  of  alignment  difficulties,  trial  runs  were  made 
with  the  tube  loosely  supported  on  a  wooden  frame,  with  the  end 
floating  on  the  vibrating  piston.  Using  this  method  of  support, 
it  was  found  that  the  tube  became  almost  self-aligning,  and  no 
other  difficulties  were  immediately  presented.  Therefore  the 
rigid  support  was  discarded,  and  the  method  of  non-rigid  sup¬ 
port  was  used  during  the  remainder  of  the  investigation. 

It  was  also  found  desirable  to  make  some  changes  in  the 
piston  as  the  experiment  progressed.  As  originally  constructed, 
the  piston  was  smooth-walled  and  made  metal-to-metal  contact 
with  the  walls  of  the  tube.  Provided  that  the  piston  was  kept 
well-lubricated  (a  mixture  of  Molykote  lubricant  and  instrument 
oil  was  used),  this  construction  seemed  to  be  fairly  satisfac¬ 
tory,  and  was  employed  for  several  weeks.  Using  this  piston, 
it  was  noted  that  changes  in  alignment  occasionally  caused  un¬ 
due  vibration  and  had  some  effect  on  the  experimental  results. 

On  one  occasion  a  combination  of  misalignment  and  insufficient 
lubrication  led  to  extensive  galling  of  the  piston  and  tube, 
and  it  was  realized  that  a  redesigned  piston  was  necessary. 

The  damaged  piston  was  machined  down  a  few  thousandths  of 
an  inch  to  remove  the  galling,  and  incidently  reducing  its  dia¬ 
meter  so  that  metal-to-metal  contact  would  be  avoided.  A  groove 
was  cut  into  which  an  O-rlng  could  be  fitted,  and  the  piston  was 
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reinstalled  on  the  shaker.  The  reconstructed  piston  proved  to 
be  more  satisfactory  than  the  old  one,  since  it  was  even  easier 
to  align,  and  was  more  nearly  gas-tight  than  before.  The  0- 
ring  was  well  lubricated  with  vacuum  grease  to  keep  wear  to  a 
minimum. 

The  end  of  the  tube  opposite  from  the  piston  was  termina¬ 
ted  with  a  movable  aluminum  plug  which  was  bored  so  that  the 
microphone  could  be  mounted  flush  with  its  face.  Two  0-rings 
provided  an  effective  seal,  and  a  long  rod  threaded  into  the 
rear  of  the  plug  allowed  the  length  of  the  tube  to  be  varied. 

d.  Oscilloscope 

A  Tektronix  Model  545  oscilloscope  was  used  as  a  recor¬ 
ding  device  because  it  possesses  certain  highly  desirable 
characteristics.  Since  this  study  involved  the  measurement  of 
steeply-rising  pressure  gradients,  a  fast  oscilloscope  rise 
time  was  essential;  the  Tektronix  scope  has  a  0.01  microsecond 
rise  time.  The  use  of  the  plug-in  preamplifier  unit  type  OA, 
which  provides  a  dual  trace,  allowed  simultaneous  observation 
of  the  pressure  signal  from  the  microphone  and  the  shaker  ve¬ 
locity  signal.  The  delay  feature  of  the  Tektronix  together 
With  its  sweep  magnifier  allowed  measurement  of  phase  differ¬ 
ences  accurate  to  about  two  degrees. 

e.  Frequency  counter 

Because  of  the  extreme  sensitivity  of  the  experiment  to 
frequency,  it  was  found  necessary  to  measure  frequency  to  the 
nearest  tenth  of  a  cycle  per  second.  For  this  purpose,  a  Hew¬ 
lett-Packard  Frequency  Counter  Model  5210  was  used.  Using  a  ten- 
second  count,  the  required  precision  in  frequency  was  obtained. 
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Figure  5-  Arrangement  of  Equipment. 


5.  Experiment  and  Results 

a.  Preliminary  Observations 

Once  the  equipment  was  set  up,  a  series  of  preliminary 
observations  over  a  wide  range  of  amplitude  and  frequency  were 
made.  The  following  are  among  the  facts  immediately  noticed: 

The  transition  from  the  acoustic  domain  to  the  shock  do¬ 
main  as  amplitude  is  increased  is  a  gradual  one.  Once  the  fre- 
quency-tubelength  combination  is  properly  adjusted  for  the  oc¬ 
currence  of  shocks,  increasing  piston  motion  from  zero  ampli¬ 
tude  produces  a  sinusoidal  pressure  output  which  gradually  dis¬ 
torts.  The  first  indication  of  an  incipient  shock  wave  is  the 
occurrence  of  a  discontinuity  in  the  first  derivative  of  the 
pressure-time  curve.  Further  increase  in  amplitude  increases 
the  bend  until  it  becomes  vertical.  Still  further  increase 
brings  about  a  break  in  the  curve  which  is  the  characteristic 
appearance  of  the  wave  containing  the  shock. 

Although  the  shape  of  the  pressure-time  wave  changes  dras¬ 
tically  during  the  transition  from  the  infinitesimal-  to  the 
finite-amplitude  region,  wave  shapes  of  almost  identical  ap¬ 
pearance  can  be  produced  at  each  frequency  for  a  given  A -value 
by  selecting  the  proper  piston  amplitude.  In  addition,  when 
the  shaker  is  oscillating  at  its  maximum  available  amplitude, 
the  pressure  wave  shapes  obtained  at  different  resonant  frequen 
cies  are  approximately  the  same.  Using  the  Altec  microphone, 
the  above-mentioned  effect  is  frequency-limited,  because  as 
frequency  is  Increased,  the  rise  time  of  the  Instrument  assumes 
a  more  and  more  significant  fraction  of  the  period  of  vibration 
masking  whatever  similarity  may  exist. 
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The  value  of  A  Is  a  very  sensitive  function  of  both  f  and 
L.  For  a  system  adjusted  so  that  A  =  0»  at  a  frequency  of, 
say,  200  cycles,  a  change  of  two- tenths  of  a  cycle  will  oause 
A  to  change  by  about  ten  degrees,  and  a  change  of  three  cycles 
per  second  will  drive  A  beyond  90°,  whereupon  the  shock  disap¬ 
pears.  Small  changes  in  the  length  of  the  tube  have  a  corres¬ 
ponding  effect  on  A  . 

Shocks  are  found  to  exist  in  the  tube  under  the  most  fa¬ 
vorable  conditions  for  their  production  (i.  e.,  large  ampli¬ 
tude  driver  motion)  through  a  A -range  of  from  about  —110°  to 
+90° .  As  A  passes  outside  of  this  range,  the  shock  diminishes 
in  size  and  finally  disappears,  although  the  vibration  does  not 
immediately  become  sinusoidal.  0 

b.  Analysis  of  Wave  Shape. 

The  fact  that  similar  wave  shapes  can  be  produced  at  each 
frequency  by  selecting  the  proper  amplitude  led  to  the  supposi¬ 
tion  that  there  may  exist  a  parameter,  independent  of  frequen- 
cy,  by  which  wave  shape  may  be  specified.  In  order  to  pursue 
this  supposition,  a  harmonic  analysis  of  the  wave  shape  was 
conducted  for  several  values  of  amplitude  and  frequency.  This 
analysis  was  made  using  a  Hewlett-Packard  Wave  Analyzer  Model 
302A.  Several  sets  of  data  were  taken  for  frequencies  of  from 
100  to  500  cps  and  for  the  second  through  the  twelfth  harmonic, 
as  functions  of  rms  piston  velocity  which  was  read  in  milli¬ 
volts  output  of  the  velocity  pickup.  These  data  were  then 
plotted  against  piston  displacement,  velocity  and  acceleration. 
It  was  soon  noted  that,  although  the  curves  of  percent  harmonic 
versus  displacement  and  versus  velocity  bore  little  similarity 
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to  each  other,  the  curves  of  percent  harmonic  versus  piston 
acceleration  were  strikingly  similar,  and  in  some  cases  iden¬ 
tical,  for  all  of  the  frequencies  tested.  It  could  then  be 
postulated  that  the  similarity  parameter  connecting  the  shocks 
is  piston  acceleration. 

In  an  effort  to  substantiate  this  postulate,  a  series  of 
twenty-seven  oscilloscope  photographs  were  taken  of  the  pres¬ 
sure-time  wave  as  seen  at  all  combinations  of  three  different 
frequencies,  three  different  values  of  A,‘  and  three  different 
acceleration  amplitudes.  In  order  to  make  a  visual  comparison, 
the  oscilloscope  vertical  sensitivity  and  sweep  time  were  ad¬ 
justed  so  that  the  wave  was  exactly  the  same  height  and  the 
same  length  in  all  cases. 

Extraction  of  the  data  from  the  oscilloscope  photographs 
was  accomplished  using  a  projection  technique.  Each  photo 
was  placed  in  an  opaque  projector  and  the  curve  projected  onto 
a  large  sheet  of  graph  paper.  Using  arbitrary  numerical  scales 
along  the  vertical  and  horizontal  axes,  the  ordinates  of  the 
pressure-time  curves  were  read  and  recorded  for  a  number  of 
equlspaced  points  along  the  abscissa.  Initially,  seventy-two 
points  were  read  off  each  curve,  but  this  was  later  reduced  to 
fifty. 

The  harmonic  analysis  itself  was  done  using  the  CDC  Model 
1604  digital  computer  installed  at  the  U.  S.  Naval  Postgraduate 
School.  A  program  was  written  which  would  accept  input  data 
in  the  form  of  the  data  points  from  the  curves,  perform  the 
Fourier  analysis,  and  then  divide  each  amplitude  by  the  ampli¬ 
tude  of  the  fundamental.  In  this  way,  the  results  were  in 
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effect  normalized  to  a  common  reference,  in  the  same  way  as 
for  the  work  done  with  the  electronic  wave  analyzer. 

The  data  obtained  in  this  way  were  then  plotted,  with  per¬ 
cent  of  the  fundamental  amplitude  as  ordinate,  and  root  mean 
square  acceleration,  A  as  abscissa.  If  acceleration  is  indeed, 
as  postulated,  a  similarity  parameter,  then  for  a  given  value 
of  A,  the  points  taken  for  common  should  be  the  same  re¬ 
gardless  of  frequency.  The  results  obtained  in  this  portion 
of  the  experiment  are  plotted  in  Pigs.  (6),  (7),  and  (8),  and 
show  that  agreement  was  obtained  to  within  at  least  two  per¬ 
cent  in  every  case,  and  thus  tend  to  confirm  the  basic  postu¬ 
late.  It  Is  believed  that  departures  from  the  anticipated  agree¬ 
ment  are  due  in  part  to  small  fluctuations  in  the  value  of  A 
for  the  points  of  different  frequency  and  to  small  changes  in 
microphone  sensitivity  to  the  different  harmonic  frequencies. 

Even  though  perfect  correlation  was  not  obtained,  it  must  be 
realized  that  if  the  harmonic  amplitudes  were  plotted  with  re¬ 
spect  to  a  parameter  other  than  acceleration  (such  as  velocity 
amplitude  or  displacement  amplitude),  the  curves  would  be  so 
dissimilar  that  they  could  not  be  plotted  on  the  same  piece  of 
paper. 

In  order  to  demonstrate  further  the  meaning  and  effect  of 
acceleration  as  a  similarity  parameter,  two  sets  of  photographs 
taken  at  different  piston  frequencies  are  included  herein  for 
comparison.  The  first  (Pig.  9)  shows  a  pair  of  oscillograms 
of  the  pressure-time  wave  taken  at  the  same  value  of  piston 
velocity  amplitude,  and  the  second  (Pig.  10)  shows  a  similar 
pair  taken  at  the  same  value  of  piston  acceleration  amplitude. 
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The  numbers  at  the  right  Indicate  the  order  of  harmonic. 
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Figure  ?.  Harmonic  Content  of 
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Figure  8.  Harmonic  Content  of 
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V  =  6.9  cm/sec  RMS 


1=6. 9  cm/sec  RMS 


Figure  9.  Oscillograms  of  pressure-time  waves  at 

different  frequencies  with  the  same  piston 
velocity. 


f  =  300  cps  f  =  500  cps 

A  -  140  m/sec 2  RMS  A  -  1  40  m/sec2  RMS 


Figure  10.  Oscillograms  of  pressure-time  waves  at 

different  frequencies  with  the  same  piston 
acceleration. 


The  reader  will  note  the  remarkable  similarity  in  the  shape 
in  the  second  case.  The  slight  difference  in  height  of  the 
shockfront  in  the  latter  set  may  be  due  to  the  fact  that  micro¬ 
phone  rise  time  assumes  a  relatively  larger  portion  of  the 
cycle  in  the  higher  frequency  case.  The  small  fluctuation 
superimposed  on  the  basic  waveform  is  due  to  the  ringing  of  the 
microphone  at  the  frequency  of  its  mechanical  resonance. 

It  must  be  pointed  out  that  merely  determining  the  har¬ 
monic  amplitudes  of  the  waves  does  not  in  Itself  give  a  com¬ 
plete  picture  of  their  construction,  since  the  complete  Fourier 
analysis  is  given  by 

P  =  A0  +  Ajcos(cot  +  ^)  +  A2cos(2ct>t  +  ♦  •  • 

and  no  determination  of  the  constants  ^  was  made  in  this  ex¬ 
periment.  Indeed,  the  fact  that  the  visual  appearance  of  the 
wave  changes  drastically  with  A  while  the  curves  of  Figs.  6 
through  8  change  only  slightly  at  larger  amplitudes  indicates 
that  there  must  be  a  substantial  change  in  the  ^  as  .A  is 
changed.  The  difficulty  in  determining  the  ^  due  "t0  the 
asymmetry  of  the  waves  which  prevents  the  selection  of  a  satis¬ 
factory  zero  reference  which  will  be  common  to  all  waves,  re¬ 
gardless  of  A •  This  same  difficulty  made  it  impossible  to  de¬ 
termine  reasonable  results  for  AQ;  the  computer  program  as 
written  solved  for  this  quantity,  but  since  no  common  vertical 
zero  reference  was  available  for  all  waves  analyzed,  the  numbers 
obtained  were  meaningless. 

c.  Variation  of  A  with  acceleration. 

In  the  course  of  the  investigation  it  was  noticed  that 
during  a  series  of  measurements  which  were  intended  to  be  at 


28 


constant  A  »  a  change  in  the  value  of At  was  often  accompanied 
by  a  change  in  the  value  of  A  ,  even  though  the  frequency  was 
maintained  constant.  Therefore,  it  was  decided  to  investigate 
this  variation.  A  series  of  measurements  was  made  which  began 
with  an  initial  setting  of  A  at  the  maximum  value  of  acceler¬ 
ation  available  on  the  shaker;  A  was  then  recorded  as  the  ac¬ 
celeration  was  reduced.  This  procedure  was  carried  out  for 
seven  different  initial  A  -settings  for  each  of  three  differ¬ 
ent  frequencies.  The  value  of  A  measured  at  the  maximum  avail¬ 
able  value  of  A  will  be  termed  A  max-  For  each  frequency  a 
similar  set  of  curves  were  obtained,  included  herein  as  Figs. 

11  through  13.  In  general,  if  the  AmQV  was  positive,  the 

uld  X 

value  of  A  became  more  positive  as  A  was  reduced;  similarly, 

a  negative  A  became  more  negative.  However,  it  was  inter- 
max 

esting  to  note  that  a  A_  of  zero  did  not  remain  zero  with 

IHd,  X 

decreasing  A  ,  but  rather  became  positive.  In  addition,  there 
was  one  negative  value  of  A  ,  which  was  different  for  each 
frequency  checked,  for  which  A  remained  constant  as  A  was  re¬ 
duced.  The  measured  positions  of  these  constant- A  lines  are 
shown  as  dashed  horizontal  lines  in  Figs.  11  through  13. 

The  close  similarity  of  these  sets  of  curves  for  the  dif¬ 
ferent  frequencies  again  demonstrates  the  existence  of  A.  as  a 
similarity  parameter. 

d.  Shock  Strength 

The  most  striking  feature  of  the  nonlinear  wave  of  finite 


amplitude  is  the  appearance  of  a  sharp  pressure  discontinuity 
as  the  shock  passes  the  recording  device.  When  the  microphone 
is  installed  in  the  rigid  tube-end  as  it  was  in  this  experiment, 


a  double-amplitude  shock  is  observed,  as  the  incident  and  re¬ 
flected  shocks  are  superimposed. 

Shock  strength  variation  was  studied  over  nine  different 
values  of  A  at  each  of  six  frequencies,  while  varying  A  over 
its  available  range.  The  strength  of  the  shock  was  measured  by 
recording  the  height  of  the  shock  as  observed  on  the  calibrated 
oscilloscope  face.  From  the  pressure  calibration  of  the  micro¬ 
phone  it  was  possible  to  obtain  shock  strength  as  a  function 
of  f,  A>  and  , A.  .  A  series  of  quasi-linear  curves  were  ob¬ 
tained  which  are  Included  herein  as  Figs.  14  through  17. 

Although  the  response  of  the  microphone  is  not  perfectly 
constant  over  all  frequencies,  it  is  sufficiently  flat  over  many 
harmonics  that  the  single  figure  of  -85  db  re  1  volt  per  dyne 
per  square  centimeter  can  be  used  to  compute  pressure  from  volt¬ 
age  output.  It  was  also  necessary,  of  course,  to  halve  the 
pressure  obtained  because  of  the  double-amplitude  nature  of 
the  wave  as  measured,  in  order  to  compute  incident  shock  over¬ 
pressure. 

The  qualitative  characteristics  of  the  shock  variation  are 
similar  for  all  frequencies  studied.  At  any  A-  setting,  as 
the  acceleration  is  reduced,  the  shock  diminishes  in  strength 
and  finally  disappears.  At  the  time  that  the  shock  disappears, 
the  wave  retains  the  form  of  a  distorted  sine  wave,  and  it  is 
not  until  is  much  further  reduced  that  the  pure  sinusoidal 
character  of  the  acoustic  wave  is  assumed.  The  value  of  A  for 
which  the  shock  disappears  is  not  the  same  for  all  frequencies, 
being  about  twice  as  great  for  1000  cycles  as  for  100  cycles. 

In  order  to  test  the  previous  conclusion  that  the  shape  of 


30 


the  wave  can  be  predicted  from  piston  acceleration,  each  value 
of  shock  strength  data  was  divided  by  the  amplitude  of  the  fun¬ 
damental  component  of  the  wave  measured  under  the  same  condi¬ 
tions.  This  was  done  to  take  into  account  the  varying  strengths 
encountered  at  the  different  frequencies.  The  points  resulting 
from  this  calculation  were  then  plotted  against  A. ,  giving  the 
set  of  curves  shown  in  Pigs.  18,  19  and  20.  These  curves  essen¬ 
tially  support  the  contention  that  A  is  a  similarity  parameter, 
since  for  every  A,  the  points  obtained  by  the  above  operation 
fall,  within  the  limits  of  experimental  precision,  along  a  sin¬ 
gle  line,  regardless  of  frequency. 

By  a  careful  study  of  the  work  of  Saenger  and  Hudson  [8j, 
one  may  conclude  that  the  existence  of  the  similarity  parameter 
A  for  shock  strength  can  be  Inferred  at  least  qualitatively 
from  their  predictions.  Their  development  results  in  the  ex¬ 
pression 
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(6) 


Pc  rr$L  (1-1) B  +  i 

where  §?  is  the  dissipation  parameter  for  viscous  losses  and 
B  is  a  parameter  derived  from  the  proportionality  constant  in 
Newton's  law  of  cooling.  According  to  Saenger  and  Hudson, 
agreement  between  experiment  and  theory  could  not  be  obtained 
unless  empirical  values  for  B  and  3?  were  used.  These  were 
such  that  3?}^B.  Hence  (6)  becomes,  approximately 
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In  a  theoretical  investigation  of  the  propagation  of  sound 
in  a  closed  tube,  Rayleigh  [9}  derived  an  expression  for  the 
tangential  force  per  unit  area  on  the  wall  by  viscous  effects 
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of  air  moving  at  unit  velocity.  Expressed  per  unit  length  of 
tube,  this  is 


F  —  TT  0.  -/  2.  Cd^o  u 


To  obtain  the  £  of  Saenger  and  Hudson  from  this,  F  must  be  di¬ 
vided  by  the  mass  of  the  gas  per  unit  length,  jOTTO.2  ( because  of 
their  use  of  the  Lagrangian  notation),  and  by  cj  (because  of 
their  use  of  nondimensionalized  quantities),  giving 


$ 


Numerically,  this  does  not  agree  with  the  empirical  results  of 
Saenger  and  Hudson.  Nevertheless,  the  important  consideration 
here  is  the  frequency  dependence;  it  is  noted  that  $  is  pro¬ 
portional  to  f“2.  Hence,  a  further  approximation  of  (7)  is 

AP  


(8) 


Applying  (5),  equation  (8)  may  be  written 

(A-#) 

which  becomes 

Sin  (A  —  $)  (9) 

Po  {Z7 JZVC  ^  v  7 

The  behavior  of  P0  is  not  well-known,  either  as  a  func¬ 
tion  of  A  or  of  frequency.  More  investigation  must  be  done  to 
reconcile  if  possible  the  facts  that,  although  AP  is  near  maxi 
mum  at  A=  0  (see  Fig.  21),  the  right  member  of  (9)  is  near 
minimum  at  the  same  value  of  A. 

e.  Phase  Diagrams 

In  taking  the  data  for  the  shock  strength  variation,  the 
minimum  value  of  J\.  for  which  a  shock  could  be  measured  was  re 
corded  at  each  of  the  resonant  frequencies  studied.  This  mini- 
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mum  occurred  at  or  near  the  A  =  0  frequency  for  the  tube. 
Plotted  against  frequency  in  Pig.  24,  these  points  furnish  a 
phase  diagram  which  separates  the  A- f  plane  into  regions  in 
which  shocks  occur  and  do  not  occur.  These  are  labelled  M Shock 
Region"  and  "Acoustic  Region,"  respectively.  A  similar  phase 
diagram,  prepared  for  displacement  amplitudes  instead  of  ac¬ 
celeration,  is  included  as  Fig.  25. 

If  A  were  in  fact  a  perfect  similarity  parameter,  then 
the  curve  of  Fig.  24  would  be  a  horizontal  straight  line. 

Since  it  is  not,  this  shows  that  although  the  concept  of  A 
as  a  similarity  parameter  for  shock  waves  is  valid  over  a  large 
range  of  piston  amplitudes,  it  breaks  down  in  the  limit  of  very 
small  amplitudes. 

The  theory  of  Saenger  and  Hudson  predicts  the  minimum 
piston  amplitudes  at  which  shocks  may  occur.  A  calculation 
based  on  this  theory  for  f^  =100  cycles  per  second  predicts  a 
minimum  displacement  amplitude  for  the  shock  region  which  is 
about  thirty  times  that  obtained  for  the  100  cycle  point  on 
Fig.  25. 

f.  Check  of  Correlation  Between  "Acoustic"  and  "Shock" 
Resonance  Frequencies. 

Although  it  has  long  been  known  that  the  production  of 
periodic  shock  waves  in  a  closed  tube  occurs  at  frequencies 
which  are  very  near  the  acoustic  resonant  frequencies  of  the 
tube,  the  writer  knows  of  no  instance  where  it  has  previously 
been  verified  that  these  frequencies  are  in  fact  the  same.  In¬ 
deed,  the  difference  in  magnitude  of  X  for  the  two  cases,  as 
well  as  the  more  complicated  theory  applicable  in  the  finite 
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amplitude  case,  leads  one  to  expect  that  they  might  conceivably 
be  slightly  different. 

A  check  on  these  frequencies  was  a  simple  matter  with  the 
variable-amplitude  driver  used.  First,  the  frequency- tubelength 
combination  was  adjusted  so  that  A  was  exactly  zero  with  the 
shaker  operating  at  maximum  acceleration.  Then  was  reduced 
using  the  GAIN  control  on  the  shaker  until  the  output  of  the 
microphone  was  near, its  minimum  detectable  value.  The  sinus¬ 
oidal  appearance  of  the  microphone  output  assured  that  this 
level  was  well  within  the  acoustic  domain.  During  the  process 
of  reducing  the  gain  control,  the  frequency  counter  indicated 
that  frequency  was  maintained  constant  to  the  nearest  tenth  of 
a  cycle  per  second. 

Finally,  the  microphone  output  was  measured  on  an  a.  c. 
voltmeter.  Several  measurements  were  made  following  this  pro¬ 
cedure,  at  frequencies  of  100,  300,  300  and  100  cycles.  In 
every  case,  the  voltmeter  was  found  to  be  sharply  peaked  when 
the  final  measurement  was  made,  and  the  slightest  change  in 
frequency  in  either  direction  caused  an  abrupt  drop  in  the 
sound  pressure  as  measured  by  the  microphone. 

From  this  it  is  concluded  that  the  acoustic  resonance  fre¬ 
quency  of  a  tube  and  the  Amax  =  0  condition  for  periodic  shock 
waves  occur  at  exactly  the  same  frequency,  within  the  limits 
of  experimental  precision. 

(In  view  of  the  shift  in  the  Amax  =  0  lines  of  Figs.  11, 
12  and  13  as  Jl.  is  decreased,  the  results  are  especially  in¬ 
teresting.  Examining  the  curves,  one  might  suspect  that  the 

values  of  Amax  for  which  =  o,  represented  by  the  dashed 

dJl 
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lines,  rather  than  the  Amax  =  ®  condition  would  be  found  to  be 
the  determining  factor  for  the  acoustic  resonance  frequencies, 
but  this  is  definitely  not  the  case.) 

g.  Variation  of  Amax  frequency. 

During  several  series  of  measurements,  the  variation  of 
Amax  with  frequency  was  studied.  As  could  be  anticipated,  it 
was  found  that  the  change  in  frequency  necessary  to  produce  a 
given  change  in  A^^.  was  greater  as  frequency  was  increased. 
The  frequency  at  which  the  tube  is  resonant  to  the  fundamental 
will  be  denoted  by  f} ;  then  the  change  necessary  to  produce  a 
given  Amax  is  expressed  by  the  quantity  f  -  f i  .  When  f  - 
was  plotted  against  A^y ,  a  family  of  monotonically  decreasing 
curves  was  obtained  (Pig.  22).  Since  all  curves  are  of  similar 
shape,  a  parameter  was  sought  by  which  the  variation  with  re¬ 
spect  to  frequency  could  be  quantitatively  determined.  After 
a  short  exploratory  process,  it  was  discovered  that  the  quan¬ 
tity  (f  -  f i  )/f]  s  is  the  same  function  of  Amax>  f°r  f  1  • 

This  result  is  shown  in  Pig.  23. 

h.  Thickness  of  Shock  Wave . 

All  oscillograms  taken  show  the  shock  wave  as  an  abrupt 
rise  in  pressure  which  takes  place  in  a  time  so  short  that  it 
could  not  be  measured  by  the  Altec  microphone  with  its  rise 
time  of  about  sixty  microseconds.  In  an  effort  to  determine 
more  about  the  structure  of  the  shock,  a  new  microphone  was 
constructed.  The  pressure-sensitive  element  of  the  new  micro¬ 
phone  was  a  barium  titanate  wafer  with  a  two-megacycle  resonant 
frequency.  The  wafer  was  mounted  in  a  brass  plug  fitted  in  the 
same  manner  as  the  plug  which  held  the  Altec  microphone. 
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Using  a  high-gain  preamplifier  with  the  Tektronix  oscil¬ 
loscope,  the  output  of  the  barium  titanate  crystal  in  response 
to  periodic  shock  waves  could  be  observed.  Although  there  was 
a  considerable  amount  of  ringing  (due  to  the  mechanical  reso¬ 
nance  of  the  crystal)  superimposed  on  the  familiar  profile  of 
the  shock  wave,  the  pressure  Jump  of  the  shock  wave  was  clear¬ 
ly  visible,  and  its  rise  time  could  be  measured.  Careful  mea¬ 
surements  at  several  frequencies  using  the  expanded  scale  of 
the  oscilloscope  gave  rise  times  of  the  order  of  one-half  mic¬ 
rosecond.  From  this  it  is  concluded  that  the  times  measured 
were  those  of  the  instrument  and  not  of  the  shock  wave,  and 
that  the  shock  pressure  Jump  must  be  completed  in  less  than 
one-half  microsecond.  The  shock  thickness,  then,  must  be 
less  than  0.2  millimeter,  which  is  about  the  distance  a  sound 
wave  travels  in  one-half  microsecond. 
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Figure  11.  Variation  of  A  as  <A.  is  reduced.  f1  =  100  cps. 

Data  was  taken  by  setting  initial  values  of 
and  recording  A  as  </£  was  reduced. 
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Figure  12.  Variation  of  A  as  JL  is  reduced,  f^  =  200  cps. 

Data  was  taken  by  setting  initial  values  of  A 
and  recording  A  as  ^2  was  reduced. 
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Figure  13.  Variation  of  A  as  <A  is  reduced.  f1  =  300  cps. 

Data  was  taken  by  setting  initial  values  of  Amax 
and  recording  as  ^  was  reduced. 


39 


40  60  80  100  120  140 

r  m  s  acceleration,  meters/sec2 


Figure  14.  Shock  Strength  versus  A. 
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Figure  15.  Shock  Strength  versus 
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Figure  17.  Shock  Strength  versus 


A. 
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Figure  18.  Normalized  Shock  Strength  versus^. 

Vertical  scale  is  Z\P  in  volts  (as  measured  on 
oscilloscope)  divided  by  fundamental  of  wave  in  rms  volts. 
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Figure  19.  Normalized  Shock  Strength  versus^. 
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Vertical  scale  is  AV  in  volts  (as  measured  on 
oscilloscope)  divided  by  fundamental  of  wave  in  rms  volts. 
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Figure  20. 


Normalized 


Shock  Strength  versus 


Vertical  scale  is  AP  In  volts  (as  measured  on 
oscilloscope)  divided  by  fundamental  of  wave  in  rms  volts. 
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Figure  22.  Variation  of  (f  -  fi )  with  ^ 
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Figure  23.  f  as  a  function  of  A 
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Figure  24.  Phase  Diagram  Showing  Threshold  Piston 
RMS  Acceleration  for  Production  of 
Periodic  Shocks 
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6.  Conclusions 

The  existence  of  acceleration  as  a  similarity  parameter 
has  been  definitely  demonstrated.  This  fact  is  important  when 
considering  the  use  of  an  electrodynamic  shaker  as  a  driver  for 
generating  shock  waves.  Such  shakers  are  rated  at  a  certain 
maximum  force  output,  which  implies  a  maximum  acceleration. 
Since  acceleration  is  a  similarity  parameter,  this  implies  that 
if  a  shaker  is  satisfactory  for  generating  shock  waves  at  one 
frequency,  it  can  be  used  for  this  purpose  at  other  frequencies 
within  its  available  range.  Because  of  the  frequency  limita¬ 
tion  of  the  Altec  microphone,  the  upper  limit  at  which  shocks 
can  be  generated  was  not  found. 

An  examination  of  the  plots  of  harmonic  analysis,  A  ,  and 
normalized  shock  strength  with  respect  to  M.  shows  a  tendency 
for  these  curves  to  level  off  at  the  higher  accelerations. 

This  fact  can  be  used  to  indicate  the  limits  of  the  threshold 
region  for  periodic  shock  waves.  Unless  unforseen  excursions 
of  these  curves  take  place  at  higher  values  of  Jl ,  it  can  be 
concluded  that  at  accelerations  of  approximately  twice  that 
available  with  the  equipment  used  in  this  experiment,  the 
curves  of  Figs.  6,  7,  8,  11,  12,  13,  18,  19  and  20  will  have 
flattened  out.  Since  the  Calidyne  Shaker  Model  A83  has  avail¬ 
able  accelerations  of  the  order  of  20g,  this  means  that  a  ma¬ 
chine  with  rated  accelerations  of  the  order  of  40g  -  50g  will 
be  of  sufficient  capacity  for  the  study  of  threshold  phenomena. 
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7.  Suggestions  for  Further  Experiment 

Saenger  and  Hudson  [8]  state  that  a  small  amount  of  ex¬ 
perimental  data  shows  qualitative  agreement  with  their  theory. 
Because  the  behavior  of  P0  was  not  measured  in  the  present  ex¬ 
periment,  no  additional  evaluation  of  the  correctness  of  the 
theory  can  be  made  here.  A  further  step  in  the  present  process 
would  be  to  obtain  the  variation  of  P0  with  A  and  A,  along 
with  shock  strength  information,  and  use  this  information  to 
obtain  a  quantitative  evaluation  of  their  theory. 

In  this  process  it  would  be  necessary  to  determine  ex¬ 
perimental  values  for  the  constants  B  and  $ .  It  is  known 
that  the  experimental  and  theoretical  values  for  B  are  in  fair¬ 
ly  close  agreement,  and  that  this  is  not  true  for  $  .  One 
possibility  would  be  to  use  the  digital  computer  to  compare 
experimental  data  with  that  calculated  using  various  values  of 
$ ,  until  a  of  optimum  size  to  suit  experimental  fact  could 
be  found.  Once  this  optimum  5  is  found,  it  could  be  deter¬ 
mined  whether  the  theory  adequately  describes  periodic  shock 
waves . 

The  relatively  slow  rise  time  of  the  microphone  employed 
in  the  present  experiment  limited  the  extent  to  which  frequency 
variation  could  be  investigated.  The  microphone  rise  time  be¬ 
came  a  significant  fraction  of  the  period  when  frequency  was 
Increased  to  1000  cycles  per  second,  and  no  Investigations 
above  that  frequency  were  considered  worthwhile.  It  is  recom¬ 
mended  that  an  investigation  be  made  of  the  effects  encountered 
at  higher  frequencies,  using  a  more  suitable  means  of  pressure 
measurement.  If  a  faster-acting  microphone  is  used  to  measure 
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pressure,  it  will  be  necessary  to  use  a  suitable  amplifier, 
with  a  filter  to  remove  the  mechanical  ringing  inherent  in  such 
equipment . 

Schlieren-optical  techniques  offer  another  means  of  inves¬ 
tigating  shock  phenomena  jjo"],  Precise  measurements  of 

density  variations  in  the  air  column  can  be  made,  and,  pro¬ 
vided  that  one  additional  parameter,  such  as  particle  velocity, 
is  measured,  the  corresponding  pressures  can  be  obtained  using 
the  equation  of  state. 
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APPENDIX  I 


SPECIFICilTI ONS  OP  MODEL  A88  SHAKER 


Characteristic 
Rated  Force  Output 
Field  Current  Required 
Continuous  Armature  Current 
Force/ Current 
Longitudinal  Resonance 
Armature  Weight 
Stroke 

Velocity  Signal  Calibration 


100  lbs  @  5.26  amp  rms 
5.0  amperes 
5.3  amperes 
19  lbs/ampere  rms 
2560  cps 
2.68  lbs 

1  inch  peak-to-peak 
0.110  volt/inch  per  sec. 
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decibels  re  1  volt/dyne  per  cur 


APPENDIX  II 


Microphone  Calibration  Curve 
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APPENDIX  III 
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TABLE  OF  SYMBOLS 

Root  mean  square  piston  acceleration 
Radius  of  tube 

Dimensionless  form  of  proportionality  constant 
In  Newton's  law  of  cooling 

Velocity  of  sound 

Driving  frequency  of  piston 

Fundamental  resonance  frequency  of  tube 

Tube  length  measured  from  mean  position  of 
piston  to  closed  end  of  tube 

Shock  strength  —  the  positive  jump  In  pressure 
across  the  shock. 

Instantaneous  pressure  at  a  point  in  the  tube 
where  gas  Is  at  rest 

Gas  particle  velocity 

Piston  velocity  amplitude 

Coordinate  of  position  along  length  of  tube 
Piston  displacement  amplitude 
Ratio  of  specific  heats,  cp/cv 

Phase  lag  of  the  piston  in  Its  outstroke  (maxi¬ 
mum  withdrawn)  position,  behind  the  appearance 
of  the  shock  at  the  piston. 

The  value  of  A  measured  with  the  shaker  opera¬ 
ting  at  Its  maximum  available  jQ.  . 

Wavelength 
Viscosity  of  gas 
Kinematic  viscosity 
Density  of  gas 

Dimensionless  form  of  viscosity  parameter 

-  [(y-i)b2-2('y-()8S  -52]/L(^-  Ob  +■  5j 
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